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Recent theoretical work on sound sources in subsonic turbomachin- 
ery is applied to the correlation of axial-flow fan and compressor 
noise measurements. Correlations are presented for directivity and 
sound-power spectra associated with broadband noise. A simple ex- 
planation is given to account for the relative amounts of sound power 
transmitted upstream and downstream from a single-stage fan. 


The mechanisms of sound generation in subsonic turbomachinery are 
understood in principle and have been documented in a recent paper by 
the author (ref. 1). Current fundamental research is now aimed at 
describing the unsteady flow within the machine with sufficient accuracy 
that the relative importance of the various mechanisms can be assessed. 

Existing noise measurements are unfortunately not matched by the 
necessary unsteady aerodynamic measurements; in this situation, positive 
identification of mechanisms is clearly not possible. Any correlation of 
data with mean-flow measurements will remain of doubtful generality 
until the gap is closed. 

If a particular noise-generation mechanism is put forward as a hypoth- 
esis, however, theory can provide considerable guidance in the correlation 
of measurements. Even though only mean-flow aerodynamic measure- 
ments are available, the relevant unsteady-flow components can be related 
to the steady flow by similarity arguments. Considerable success has been 
achieved in this way in jet noise correlation, following the initial work of 
Lighthill (ref. 2) ; the same approach is adopted in this paper for broad- 
band fan noise. 


1 The cooperation of Rolls-Royce Limited in supplying data is gratefully acknowl- 
edged. The initial calculations were carried out by R. Snow, Il.A.C. Acoustics Labora- 
tory, Weybridge, while at I.S.V.R., Southampton University. 
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Following the studies by Smith and House (ref. 3) and Morfey and 
Dawson (ref. 4) on aircraft-engine fan and compressor noise, further data 
on broadband noise have become available covering a range of axial-flow 
machines (table I). These are reported by Snow in reference 5, which 
gives details of measurement and analysis techniques and forms the basis 
of the present study. 

In the sections that follow, various aspects of the noise measurements — 
sound-power spectrum shapes and levels, directivity patterns and the 
ratio of downstream to upstream radiated power — are discussed and 
correlated with mean-flow data. The theory on which the correlations are 
based has been developed where necessary in the appendixes. 


CORRELATION OF BROADBAND POWER SPECTRUM SHAPES 

The correlation of broadband sound-power spectra is approached in two 
stages. First, a frequency parameter is sought that brings tin' different 
spectra into lint' horizontally; then, for a fixed value of the frequency 
parameter, normalized values of sound-power spectral density are plotted 
versus Mach number. This is the same general procedure as was followed 
by Smith and House (ref. 3), although the details differ. 


Noiiclimeusioual Frequency Parameter 

Given a characteristic length scale L, a nondimensional frequency 
parameter can be formed either as a Strouhal number fL/U, where U is 
a characteristic velocity, or as the ratio fL/c ( = L/\). These are related 
by the Mach number U/c. 

We start by choosing a length scale, leaving till later the choice between 
V and c as velocity scales. It is expected that the spectrum of broadband 
sound radiated from a rotor-stator stage will depend on the scale of 
turbulence in the rotor wake. This in turn is related to the rotor wake- 
momentum thickness, as Mugridge has shown (ref. 6). 

A possible length scab' is therefore b R Co (equal to twice the momentum 
thickness), where b K is the first-rotor chord and Co is the profile drag 
coefficient. This may conveniently be estimated in terms of the irreversible 
temperature rise A T irr through the first rotor (ref. 7) : 


b R C D = 


(A Th 

To, 


2 cIr cos 02 

(V-1)M K * 


(l ) 2 


In equation (1) , d R cos /?2 represents the rotor wake spacing and Mr is 
the relative Mach number at the rotor exit. 


2 (A7’) irr is assumed small compared with the stagnation temperature T 0 , at exit 
from the first rotor. 
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For the engines tested, information was available on the overall 
isentropic efficiency 77 ,, and hence on the irreversible temperature rise 
(1 — 77 ,,) (AT 0 )„ over all n stages. Neglecting stator losses and appor- 
tioning the irreversible temperature rise equally between the n rotors 
gives as a rough estimate 

{AT) itr = l (l~ Vn )(AT 0 ) n (2) 

n 


Combining equations (1) and ( 2 ) gives a length scale L, defined as 


2(1— 7?,,) (AT'o)n 

n(y-l)T 01 M R 2 


•d K cos ft 


(3) 


which because of its link with the wake-momentum thickness should 
reveal any variations in the eddy scale leaving the first rotor. 

This length scale is used below to obtain an empirical collapse of broad- 
band spectrum shapes, with one modification: The rotor relative Mach 
number is taken for convenience as the inlet (rather than the exit) value. 
Midradius values are used throughout. 


Collapse of Spectrum Shapes 

Figure 1 shows, for each engine separately, the collapse of inlet sound- 
power spectra obtained with L/X = F as frequency parameter. Somewhat 
surprisingly, this is a good deal better than the collapse obtained with the 
Strouhal-type parameter F/M R . Note that the spectra are plotted in 
constant-percentage-bandwidth form (i.e., the vertical scale is propor- 
tional to f(dW/df ) ) ; the peak occurs around F = 0.1 to 0.2 in each case. 
No correction has been made for atmospheric attenuation; the greatest 
measuring distance was 32 feet. 

In figure 2 the mean spectrum shapes from the different engines are 
superimposed; the broken line for engine B indicates that the horizontal 
positioning of the spectrum is uncertain since the efficiency had to be 
estimated from other data. For comparison, figure 2 also shows the mean 
spectrum shape obtained by Smith and House (ref. 3) for broadband com- 
pressor noise at the position of maximum sideline OASPL; the peak has 
been placed at F = 0.16 to line up with the present data . 3 

A mean line based on figure 2 lies within 5 dB of most of the data points. 
Some of the scatter is doubtless due to errors both in sound-power meas- 
urement and in the aerodynamic data used for forming the frequency 
parameter; but the present scheme cannot account for genuine variations 
in spectrum shape. 


3 This implies a value of Co between 0.1 and 0.2 for the compressors plotted in 
figure 12 of reference 3 
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Figure 1. — Spectrum shape correlation — broadband sound power radiated from inlet. 
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F 

Figure 2. — Mean spectrum shapes for all engines ( constant-percentage bandwidth). 


CORRELATION OF ABSOLUTE LEVELS 

The next stage in the data correlation process is to collapse the sound- 
power spectrum levels for a given value of the frequency parameter. This 
will be done with the aid of a simplified analytical model that predicts the 
sound-power variation with Mach number, blade area, and other prin- 
cipal factors. A nondimensional sound-power parameter is obtained in 
which these factors are normalized out; it is found to be constant within 
±5 dB for all the engines measured, up to a rotor inlet Mach number 
of 0.6. 

Nontlimensioual Sound-Power Parameter 

The broadband sound-power spectrum generated by a single rotor- 
stator stage is related in Appendix I to the characteristics of the rotor 
wake turbulence. The analysis is highly simplified, its purpose being to 
predict the main trends so that these can be normalized out from the data. 
The analysis is based on the following assumptions: 

(1) Broadband noise is generated by turbulence from one blade row, 
entering the row immediately downstream and setting up lift fluctuations 
on the downstream blades. Only the axial component of the lift fluc- 
tuations is taken into account in calculating the sound output (this point 
will be mentioned again below) . 

(2) Blade profiles are assumed acoustically compact. (This is a rather 
crude assumption since the ratio of axial chord to wavelength is of order 
1 at the peak of the spectra in figure 1. The effect of noncompact blade 
profiles on sound radiation requires investigation.) 

(3) The axisymmctric rotor-stator geometry is replaced by a two- 
dimensional model based on mean-radius conditions, and the spanwise 
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extent of the model is taken as being large compared with the sound 
wavelength. 

(4) The flow is subsonic with M 2 « 1. 

The results of the analysis are summarized in equations (15) and (16) 
of Appendix I. These involve the unknown quantities l c (spanwisc integral 
length scale of turbulence) and C D (rotor-blade drag coefficient) , besides 
the turbulence intensity parameter a. However, if we assume l c « b R C D , as 
Mugridge’s experiments (ref. 6) indicate, and replace b R C D by L as in the 
previous section, we get the following relationship between the sound- 
power spectrum and aerodynamic parameters : 

t) ° snM ' a ( F ) w 

The left-hand side of equation (4) is itself a nondimensional form of 
the sound-power spectrum; it will be divided by the first four factors on 
the right-hand side to produce a normalized power parameter. 

For this preliminary study, it was considered an advantage to base the 
normalization solely on rotor parameters; <x s is therefore replaced by the 
rotor solidity a R , and a\ is replaced by ft. Furthermore the C L term in 
/ (ft,«i) is neglected. The resulting nondimensional sound-power parameter 
is 


dW Pd R * es c 4 ft 

df ' p&Sn'Mf? 1 

Here Sr' = <trSr denotes the rotor blade area (span X chord X number of 
blades). All the quantities that make up G, apart from the sound-power 
spectrum itself, are generally known at the design stage of a fan or com- 
pressor. 

As in the previous section, the correlation that follows is based on the 
relative Mach number at the rotor inlet. 4 

Sound Power Data Correlation 

Values of G for F = 0.16 — corresponding approximately to the peak of 
the constant-percentage-bandwidth spectra in figure 1 — are plotted versus 
the rotor inlet relative Mach number in figure 3. The data points represent 
inlet-radiated broadband sound from the six different engines in table I, 
each run over a range of speeds. 


4 Note that for a given F value, G cc Mr *; so the value of G is not very sensitive to 
the exact definition of Mr, provided the same definition is used in both F and G. 
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Figure 3. — Sound-power parameter based on peak of spectrum, versus first-ro'.or inlet 

relative Mach number. 


Ideally the sound power from each stage would be correlated separately 
when dealing with a multistage machine. However, there is no way in 
which the broadband sound output from separate stages can be dis- 
tinguished when only the total sound output has been measured, as has 
been done here. The inlet sound-power measurements from the multistage 
machines (A, B, C, D) have therefore been combined with first-stage 
rotor parameters to give a value of (!. Such a normalization procedure 
leaves the number of stages unaccounted for, so it is interesting that the 
single-stage fan tends to fall at the bottom of the scatter band in figure 3 
(i.e., about 5 dB below most of the multistage points). 

Below = 0.6, almost all the points lie within ±5 dB of the constant 
value 6' = 0.03; this collapse is encouraging in view of the small number of 
variables included in the definition of G. Furthermore, the tendency to 
remain constant over this range of Mach number supports the theoretically 
predicted A/ft 5 variation (see equation (4)). 5 


6 Note that if the frequency spectrum scaled on relative velocity, rather than the 
sound speed as indicated in figure 1. the right-hand side of equation (4) for l c «\ would 
he proportional to 4 /k 6 rather than M R h . Although Strouhal number scaling— and 
hence an Mr* dependence for compact sources — is commonly assumed for sound radia- 
tion from fluctuating forces of aerodynamic origin, we have shown (“Correlation of 
Broadband Power Spectrum Shapes”) that it does not fit the present data. 
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Above Mr = 0.G, figure 3 shows a falloff in G with increasing M R \ this 
amounts to roughly 10 dB at M R = 0.7 and is similar to the flow effect 
found by Smith and House (ref. 3) . 

DIRECTIVITY OF INLET-RADIATED SOUND 

The directivity of broadband inlet-radiated sound is discussed in 
reference 8; the opportunity is taken here to present an updated correla- 
tion due to Snow (ref. 5). In the light of these results, a hypothesis is 
put forward to account for the directional radiation of broadband sound 
at high frequencies. 

Multistage Data Correlation 

Figures 4a and 4b show the directivity index, 6 plotted versus angle from 
the inlet axis, for five different multistage fans and compressors (n= 2 to 
8). The parameter is ka, where a is the rms first-rotor radius; each ka curve 
is based on a number of engine speeds. 

At a nondimensional frequency of kn = 4, the curves for the various 
engines collapse within 1§ dB. They also agree fairly closely (up to 60°) 
with the theoretical curve in figure 5. The latter is taken from reference 
8 and represents a spatially random source distribution; the departure 
from a uniform field shape is due to duct cutoff effects. 

At ka= 16, the same theoretical model predicts a fairly flat directivity 
curve, with a slight dip at 90° (fig. 5) ; but the measured curves fall off a 
good deal more sharply away from the axis. The departure from the 
random-source model is still more pronounced at ka = 63. In fact, the 
model predicts hemispherically uniform radiation in the high-frequency 
limit (ka = <* ) , while the measured radiation patterns become progres- 
sively more directional with increasing frequency above ka = 4 and 
apparently tend toward a cos 2 \p directional distribution. 7 

Optical Model for High-Frequency Directivity 

A clue to the observed directivity along the axis at high frequencies 
(ka= 16 and above) is provided by closer examination of figure 4. There 
is a significant trend at 90° for the directivity index to vary according to 
the number of stages: The higher the number, the larger the dip at 90°. 
An optical model, outlined below, suggests that reflection of sound by 
upstream blade rows could account for the observed directivity patterns. 


6 Bused on the upstream hemisphere (^=0° to 00°); thus omnidirectional radiation 
has a directivity index of 3 dB. 

7 It is interesting that this is the distribution which corresponds to fluctuating forces 
oriented in the axial direction. 
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Figure 4a. — Narrowband directivity index measurements for inlet radiation from multi- 
stage fans and compressors ( frequency band width: 5 percent). 


In a multistage machine, sound generated by the second and subsequent 
stages must travel through the first, stage before it can escape from the 
inlet. The transmission process is represented in figure G by rays that are 
reflected from the blade surfaces in the same way as light from a mirror; 
this optical model becomes valid at high frequencies, where the sound 
wavelength is short compared with the blade chord and spacing. 

Referring to figure 6, rays incident on the stator within the shaded band 
will be reflected back in the downstream direction. The same applies to 
the rotor, so in this two-dimensional picture the radiation transmitted 




Figure 4b. — Narrowband directivity index measurements for inlet radiation from multi- 
stage fans and compressors (frequency bandwidth: 5 percent). 


upstream is confined to a wedge of included angle tt — 2(a 2 '+ft>')- 8 Return- 
ing to the three-dimensional geometry of the actual rotor-stator stage, 
the optical model clearly predicts a reduction of transmitted intensity at 
90° (compared with the axial direction) due to upstream blade rows. 

This effect probably accounts for at least part of the observed direc- 
tivity of inlet-radiated noise at ka= 16 and above. The similarity in figure 
4 between the fans A (two-stage) and B (four-stage) may be connected 
with the absence of IGV’s in the latter. 


8 An interesting point is that, if the rotor and stator blade angles add up to more 
than 90° (which is unusual in aircraft-engine fans and compressors), the model pre- 
dicts zero sound transmission through this stage. 


E 
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Figure 5. — Theoretical directivity index based on random source model, for ka =0 to x. 
Hands indicate hob-tip ratio range 0.S-0.5. 


*2 



ROTOR ST A TOR 

Figure C. — High-frequency sound transmission through blade rows. 


EFFECTS OF SOURCE MOTION AND MEAN FLOVt 

The relative amounts of broadband sound power escaping from the 
exhaust and inlet of a single-stage' fan have been studied by Smith and 
House (ref. 3) ; they show an empirical correlation with the' relative Mach 
number through the rotor. The same measurements have been reanalyzed 
for this study tei give the downstream upstream peiwer ratio in narrow 
frequency bands. The tendency reported in reference 3 for more sound 
power to travel downstream is verified, but it appears to be restricted to 
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high frequencies (above 2 kHz). The high-frequency data are compared 
below with the predictions of a theoretical model that accounts for source 
motion and mean flow. 

Theoretical Prediction of Downstream /Upstream Energy Split 

The single blade-row model of Appendix III gives an estimate of the 
sound-power spectrum radiated downstream and upstream due to broad- 
band blade forces. Provided the blades are not highly staggered (mean 
relative flow angle less than 30°), the downstream/upstream power ratio 
can be approximated to first order in Mach number by 


1 + \M X 
1-f M x 


(6) 


The total sound-power spectrum (inlet+exhaust) is unaffected by flow 
to this approximation. 

The energy split predicted by Appendix III actually depends on the 
shape of the radiated power spectrum, and equation (G) is strictly valid 
only for a flat spectrum of sound power in constant-percentage bands. 
Furthermore, the analysis is based on the assumption that many modes 
contribute to the sound-power spectrum at the frequency concerned. 
These considerations limit the application of equation (G) to frequencies 
around the peak of the constant-percentage-bandwidth spectrum. 


Comparison of Theory with Experimental Data 

Figure 7 shows the downstream/upstream power ratio for engine E, 
measured at three different frequencies (3.8, 7, and 11.5 kHz) over a 
range of axial Mach numbers. The constant-percentage broadband spec- 
trum peaks at around 4 to 7 kHz at all speeds. 

Although the data are scattered, there is a definite trend toward in- 
creasing power ratios with increasing axial Mach number, a finding which 
agrees fairly well with equation (G) (shown on fig. 7 as a solid line). 


CONCLUSIONS 

(1) Measurements of inlet-radiated broadband power spectra from a 
range of fans and compressors show a collapse of spectrum shape based on 
the estimated first-rotor wake-momentum thickness (figs. 1 and 2). 

(2) A simplified model of broadband sound-power generation has ueen 
used to establish the main factors determining the radiated power, leading 
to the definition of a nondimensional sound-power parameter in which 
these factors are normalized out. 
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Figure 7. — Ratio of exhaust to inlet sound power in narrow frequency bands — single- 

stage fan. 

(3) Plotting the peak (constant-percentage) power spectrum levels 
in normalized form collapses the data within 5 dB (fig. 3). Below a rotor 
relative Mach number of 0.6, the nondimensional power parameter is 
approximately independent of Mach number; above 0.6 there is a notice- 
able falloff, similar to the flow effect found by Smith and House (ref. 3). 

(4) The single-stage fan levels tend to fall at the bottom of the scatter 
band; i.e., about 5 dB below most of the multistage points. This may well 
be a stage-number effect; figure 14 of reference 3 predicts differences of 
this order between single- and multi-stage machines over the relative 
Mach number range concerned (0.2 to 0.6). The difference largely dis- 
appears if the exhaust-radiated power is included in the single-stage data. 

(5) Inlet directivity measurements in narrow frequency bands show a 
consistent trend toward more directional radiation at higher frequencies 
(fig. 4), which becomes more marked the larger the number of fan or 
compressor stages. The high-frequency beaming along the axis is, if any- 
thing, more pronounced at low engine speeds (jet noise obscures the 
pattern at the highest speeds) and thus cannot be explained in terms of 
refraction. Transmission through upstream blade rows provides a qualita- 
tive explanation. 
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(6) For a single-stage fan, the ratio of exhaust to inlet sound power 
has been measured in narrow frequency bands. Around the peak of the 
constant-percentage broadband power spectrum, the trend of the measure- 
ments agrees with the ratio [1+ (fd/ I )]]/[l— (jM x ) ] theoretically 
predicted for low-stagger blade rows (fig. 7). 
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Vn 
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rms radius, outer radius of annulus 
blade chord 

blade drag, lift coefficients 
tangential spacing of blades 
frequency parameter (L/X) 
frequency 

defined in equation (9) 
sound-power parameter, equation (5) 
frequency exponent, equation (23) 
acoustic intensity 
acoustic wavenumber (w/c) 

turbulence wavenumber components (tangential, spanwise) 
length scale, equation (3) 
axial blade force per unit span 

spanwise correlation length of turbulence, equation (11) 

rotor Mach numbers (relative, tangential) 

flow Mach numbers (axial, peripheral) 

annulus mode numbers (circumferential, radial) 

number of stages in machine 

area of duct cross section 

rotor disk area, blade area 

absolute temperature, stagnation value 

overall stagnation temperature rise 

sound power, spectral density 

defined in equation (19) 

turbulence intensity parameter, equation (13); axial wave- 
number parameter, Appendix III 
stator flow angle, blade angle 
rotor relative flow angle, blade angle 
specific-heat ratio 

ratio of tangential to axial force at outer wall of duct 

overall isentropic efficiency 

sound wavelength 

fluid density 

blade-row solidity 
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\p angle from duct axis 

co radian frequency (relative to swirl in Appendix IV) 

wo, «' radian frequency in fixed, rotor reference frames 

Subscripts 

R, S rotor, stator 

1, 2 inlet, exit values for blade row 

( ) denotes time average 

APPENDIX I 

Generation of Broadband Acoustic Energy 

Broadband noise generation is here attributed to the interaction be- 
tween flow turbulence and the fan or compressor blades. More specifically, 
the process of sound generation at subsonic speeds is viewed in two stages: 
Turbulent flow produces lift fluctuations on a row of blades, and the lift 
fluctuations give rise to sound radiation. We shall consider the particular 
case where the turbulence is produced by the row of blades immediately 
upstream, although this is not the only possibility. The following analysis 
is designed to give guidelines for data correlation and is by no means a 
complete description of the complex processes involved. 

Periodic Interaction Theory 

A relatively simple case of rotor-stator interaction is that associated 
with the mean rotor-wake profile. Reference 1 gives the sound intensity, 
at any multiple' of blade-passing frequency, due to a uniform axial force 
L x per unit span on each stator blade' at that frequency: 

/ = — --- - (single mode, well above cutoff) (7) 

4 pc d s 2 

Equation (7) is based on a two-dimensional model and assumes low- 
Mach-number flow. With the same assumptiems, the mean-square har- 
monic force is given by an extension of the Kemp-Scars analysis (re'fs. 
1 and 9) as 

(Lx 2 ) = l K u . P — a,) (8) 

1 CO 

where K w is a facteir allowing for the- downstream decay of the rotor 
wakes and 

cos Bt [ Cr 

= tan 2 a\ '■ sill 2 (ft + exi) ■ 1 cot (fti + ai) 

cos L tt 





(9) 
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Note that the drag coefficient Cd refers to the rotor blades, and the lift 
coefficient Cl to the downstream stator blades. The blade-row geometry 
is shown in figure 8. 

Application to Broadband Turbulence 

The single-mode, single-frequency results given above correspond to a 
rotor-wake pattern that is sinusoidal in the tangential direction (i.e., 
single Ay wave number), uniform in the spanwise direction (Ay = 0), and 
carried round with the rotor. To deal with incident turbulence that has a 
continuous range of wave numbers and frequencies, we use the wave- 
number-frequency spectrum is (Ay, Ay, o) of the fluctuating velocity 
entering the stator. Only velocity components parallel to the rotor wakes 
(i.e., streamwise relative to the rotor) are considered here. 

Figure' 9 shows which regions of the turbulence spectrum give rise to 
propagating acoustic modes. The criterion for propagation is that for a 
given radian frequency oi the interaction wave number — whose com- 
ponents are (k y +2ws/d s , Ay), where s= dbl, ±2, etc. — should be less than 
the acoustic wave number k = u/c. The contributing regions in the (Ay, Ay) 




Figure 8. — Geometry of rotor-stator stage. 
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Figure 9. — Regions of (k y ,k z ) wave-number spectrum of turbulence which contribute to 
sound radiation at frequency u. 


plane are therefore circles of radius k ; these overlap if the stator blades 
are spaced further apart than the acoustic wavelength, and are distinct 
regions if ds<\/ 2. 

In either case, the equivalent of the mean-square harmonic wake velocity 
in the previous section is an integral of the wave-number spectrum over 
all the contributing regions. However, exact correspondence occurs only 
in the center of each circle; departure from two-dimensionality and 
proximity to cutoff introduce a weighting factor for points away from the 
center proportional to the axial wave number of the radiated sound 9 and, 
hence, proportional to the height of the hemisphere with the circle as base. 

We make the following assumptions in order to estimate the wave- 
number integral: (1) The variation of E(k y ,k z ,w) is small over an interval 
2ir/d s in fc„; (2) the variation with k z is small over the contributing region 
(| k z |</c).Thus 

E(k y ,k z , w) = E (k y> 0,u) 


This implies that the integral length scale l c (see below) is small com- 
pared with X. 

With these approximations, and taking account of the weighting factor 
mentioned above, the contributing region of the ( k y ,k z ) turbulence 
spectrum is estimated as 


2 

3 




E(k y ,0,u) dk y 


( 10 ) 


(Note factor 2 for + and — frequencies.) 


9 This is an acoustic coupling factor for axial forces (see ref. 10, sec. 6.4). The 
aerodynamic effect of finite k z on the fluctuating lift is shown in Appendix II to be 
unimportant. 
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If the integral length scale l c in the spanwise direction is defined by 


l c = 2irE (k y ,0,u) 


E ( k y , A'2 jCo) dl\y Z 


(11) 


and the dependence of l c on k y is neglected, expression ( 10) may be written 
as 


— k 2 dsl c G(u ) 

Dir 


( 12 ) 


In the above expression 


G( co) =2 


/ (» r 00 

dky / 

-00 ^ — oc 


dk zE (k y jk z j Co) 


is the single-sided co-spectral density of the velocity fluctuations entering 
the stator. A nondimensional measure of G(w) is provided by the ratio a 
of uG(u) to the mean-square fundamental 10 component of the rotor-wake 
velocity defect: 


a (a>) = ccG(co) • 


2 

cwCdW 7 2 


(13) 


It follows from equations (12) and (13) that, in order to apply the 
periodic-wake analysis to the present problem, it is simply necessary to 
replace K w in equation (8) by 

-‘-fc 2 d s lr- (14) 

07T CO 


This, with equation (7), gives the acoustic intensity per unit bandwidth 
due to axial lift forces on the stator blades as 


dj_ 

dco 



pcHcasa^C^M^f (fo, ai) 


(15) 


with ? c «\. 

The opposite extreme to equation (15) is obtained when the spanwise 
integral length scale is greater than X ; a similar analysis then gives 

~ = ape 3 • - WC'sW/^aO ( 16) 

aco 62 co 


with l c >\. 


10 i.e., with period equal to the rotor blade space. 
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Radiation from Turbulent Boundary-Layer Pressure 
Fluctuations 

In the preceding analysis, blade lift fluctuations are attributed to 
incident turbulence. An alternative cause of lift fluctuations is the un- 
steady pressure field sot up by the blade’s own turbulent boundary layer 
and wake. This mechanism was recognized by Kramer (ref. 11), and a 
quantitative estimate of t ho sound radiation was made by Sharland (ref. 
12 ) using theoretical ideas put forward by Lilloy (ref. 18). In fact the 
simplified sound power parameter, equation (.’>), used to correlate the 
data does not distinguish between the two mechanisms. The question, 
whether the rotor-generated turbulence radiates by interaction with the 
rotor or the downstream stator, is therefore left open. 


APPENDIX II 

Three-Dimensional Effects on Fluctuating Lift 

Mugridge (ref. 6 ) has estimated the fluctuating lift per unit span on a 
two-dimensional airfoil in a three-dimensional converted gust pattern by 
considering a single radian frequency u and spanwise wavenumber k,. 
The case k z = 0 corresponds to Sears’ two-dimensional analysis (ref. 9). 

The mean-square fluctuating lift is approximately equal to the two- 
dimensional mean-square lift multiplied by a correction factor 

with p 2 +n 2 < 4, where n = wb/2U (i.e. , the reduced frequency) and 
V = 2 kj>. 

Since p/n = k z U/oi, and we are interested only in values of k z less than 
k for sound radiation, it follows that over this range p < Mn with M = U/c. 
This limits the correction factor above to between 5 and 1 for subsonic 
flow. 

For high gust wavenumbers, a more accurate value of the correction 
factor is available from an independent analysis by l'ilotas (ref. 14) ; it is 

— ( 1 +P 2 n 2 )-' /2 

irp 

with p 2 >4, n >1. 

Large values of p will be associated in the present problem with acoustic 
wavelengths that are small compared with the blade chord (i.e., high fre- 
quencies). However, the acoustic representation of the blade as a line 
force breaks down in this region. 
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APPENDIX III 


Downstream/Upstream Energy Split from an Isolated Blade Row 

Apart from any reflecting surfaces upstream or downstream, a rotating 
blade row is expected to radiate different amounts of sound power in the 
upstream and downstream directions. Both the source motion and mean 
flow contribute to this effect, which applies both to discrete-frequency 11 
and broadband radiation. The following analysis deals with the broad- 
band case. 

For this purpose the blade row is modeled bv a force distribution at one 
cross section of a uniform axisymmetric duct. The local force per unit area 
has components (f x ,fo,fr) in the axial, tangential, and radial directions, 
and we assume for simplicity that 

/» = <(£)/x (17) 


with/ r = 0. (ao is the outer radius of the annular cross section.) A swirling 
mean flow in the duct is modeled by a uniform axial Mach number M x 
together with a solid-body rotation. 

The analysis starts from the result derived in reference 16, section 4, 
for single-frequency excitation in a single (m,n) duct mode. This gives 
the downstream (W+) and upstream (IF-) sound power in terms of the 
axial-force modal amplitude F x , as 


W±(m,n ) 


1 <*> | F z |*5 . 
1 — A 

8 w ape 


and 


A± = 


±a-M,+ (coo/'oi) (mt/w)(\ — Mx) 2 
lTaM x 


(18) 


Here a is the axial wavenumber parameter, defined in equation (24) 
below; wo = iec/ao is the fixed-frame radian frequency, and w is the corre- 
sponding frequency relative to the swirl. Equation (18) is valid over the 
whole frequency range above mode cutoff. 


Frequency Relationships 

If the swirl angular velocity is fi in the 6 direction, « is related to wo by 
w = w o — mSl {m = circumferential mode number) 


11 The discrete-frequency case has been studied theoretically and checked against 
experimental data by Mani (ref. 15). 
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Thus 


W „ 171 ,r 

— =1 Alg 

W 0 IV 


( 19 ) 


where 


w = — 0 (based on fixed-frame frequency) 

c 


and 



is the peripheral swirl Mach number. 

Similarly the radian frequency w' in the rotor frame of reference is 
related to w 0 by 



m 

- M t 


w 0 w 


(20) 


where M , is the rotor tip Mach number. 


Broadband Excitation 


Since w 0 and w' arc related by a constant factor for a given m value, the 
single-frequency sound-power result (eq. (18)) can be generalized to 
broadband excitation by the substitutions 


Thus 


i | F x | WG x (w') 


W±(m,n)—> wo 


dW±(m,n) 

(iwn 


dn ± l o w 

= T b*(w )‘A. 

a wo 4 ape w 


( 21 ) 


( 22 ) 


If the force spectral density G x { w') varies as (w')°> equation (22) can be 
written as 


dW ± (m,n) 

dwo 



G x (wo ) • A± 


(23) 
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Cutoff Condition 


We consider only (m, 0) modes, for which 





Cutoff occurs when 


(24) 



(1 — M* 2 ) = 1 


i.e., from equation (19), 



w 


2 


Thus the propagation range is given by 

— Mp— (1— AT, 2 ) 1/2 m -M x 2 )' /2 

1 -M 2 ~w~ 1 — M- 


(25) 


with M 2 = My + M t 2 . diodes within this range will propagate along the 
duct; modes outside the range make no contribution to the sound power 
in the present model. 

The substitution 


m 

iv 


-Me (1 
1—M 2 ' 1 -M 2 


sin <j> 


(26) 


gives the propagation range in terms of </> as — Also, a is 

given by 

a = — ( 1 — M 2 ) 1/2 ( 1 — M x 2 ) 1/2 cos 0 (27) 

CO 


Summation over Propagating Modes 

The total sound power per unit bandwidth dW ± /dun >, radiated in either 
direction, is found by summing the contributions given by equation (22) 
for all propagating modes. Two simplifying assumptions are introduced at 
this point. 

(1) The excitation spectral density G x is taken as the same for all 
(ni, 0) modes and as zero for all higher order radial modes. Although it 
would probably be more realistic to assign equal excitation to all modes 
regardless of radial order, the exclusion of higher order modes is not 
expected to alter the downstream/upstream energy split significantly. 
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(2) The propagation rango is assumed to contain a large enough 
number of (m,0) modes that the discrete summation over integer m 
values can be replaced by integration over a continuous range. The multi- 
mode sound power in either direction is then 


clW ± 

du>o 



dW±(m,0) 


daio 


dm 


w(l — M x 2 ) 1 ' 2 dll ± Q.,0) 

1 — -I -x, 2 da> o 


| 

[ 

cos <p-d<t> ) 

J 


(28) 


The integral over <f> in equation (28) is a function of e, M X) Me, and M,) 
it can be evaluated using equations (18), (23), (2G), and (27). A particu- 
larly simple result, valid for low Mach numbers, is obtained by neglecting 
squares and higher powers of the various Mach numbers; to this accuracy, 


d\V ± 
d o) 0 


- —w 6 Y x (wo) 
4 pc 


j7T (1 + e 2 ) T - [(1 — e 2 ) .V, + « .1/j + (<7 + 1) tM (29) 


Special Case: Force Perpendicular to Relative Mean Flow 

If the fluctuating blade forces act in the lift direction (i.e., perpendicular 
to the relative flow), then 

M 

e = cot dm — — — - — (30) 

M t — JV1 e 

where dm is the mean relative flow angle referred to the axis. Combining 
this relation with equation (29), and taking the exponent g as —2, gives 

f/ir+/ f/coo 1+ {8/3ir)M x COS 2 dm . 

dWI/d^o ~ W 8/3 tt)M x cos 2 dm 

The assumption g = — 2 corresponds, according to equation (29), to a flat 
spectrum of sound power in constant-percentage bands. 

A noteworthy feature of equation (31) is that, provided dm is not more 
than 30°, the downstream/upstream sound power ratio is almost in- 
dependent of dm and can be approximated by 

1 + 1 M* 

T-\m x 


( 32 ) 
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In other words, the energy split is determined solely by the axial Mach 
number and is the same for both rotating and stationary blades, provided 
the relative flow angles are not too large. 
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DISCUSSION 


M. V. LOWSON (Loughborough University, England) : This inter- 
esting paper by Dr. Morfey raises a number of points which deserve 
further discussion. It is assumed in the paper that shed turbulence is the 
principal source of broadband noise radiation. While this is certainly an 
important contributor, other sources can also give rise to broadband 
noise. The turbulence in the inlet boundary layer of a jet engine is prob- 
ably much larger, in both scale and level, than that of the standard 
equilibrium boundary layer, and also interacts with the highest speed 
part of the rotor near the tip. Indeed the tip itself can be regarded as 
running in a stalk'd condition, and the direct self-induced radiation at the 
tip could well be large. Definitive measurements of inlet turbulence 
levels and spectra are badly needed to illuminate this problem. 

The empirical correlation on momentum thickness demonstrated by 
Morfey offers some evidence in favor of the dominance of shed turbulence 
effects, but it should be pointed out that, because all engines are designed 
to similar rules, empirical parameters tend to move together. Thus there 
is some danger in inferring source mechanisms from the success of em- 
pirical collapses. 

The measured high-frequency beaming of the broadband noise is an 
interesting effect, particularly since it is more pronounced than pre- 
dicted by the theory. A possible reason for this is that tin 1 theoretical 
estimates ascribe too high a level to the higher order circumferential and 
radial modes. The actual source strength must fall off smoothly toward 
the tip and this would cause a lower level of the higher order radial modes 
than the rectangular distribution assumed. Since it is these higher order 
modes which cause the radiation normal to the axis, a somewhat lower 
level of sideline radiation than that of the present theory can be an- 
ticipated. 

As a final point, I wonder if Dr. Morfey would care to indicate his 
present thoughts on the likely relative significance of the dipole and 
quadrupole source mechanisms for broadband noise radiation. 

M. E. HOUSE (Southampton University) : Dr. Morfey is well recog- 
nized for his theoretical contribution to the problems of turbomachinery 
noise, and indeed his paper presented at the 1969 Winter Annual Meeting 
of ASME held at Los Angeles, California (ref. 1), contains a very full 
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treatment of subsonic rotor-stator interaction noise in terms of steady- 
state parameters. I must confess, however, to being a little disappointed 
that that paper offered no experimental verification. 

His present paper extends the concepts to the case of random or broad- 
band sound, and here I am pleased to note that experimental data are 
included. 

I have little or no comment to make on the derivations in the paper or 
the appendixes since they follow logically from the basis and assumptions 
stated by Dr. Morfey. My main concern lies within these assumptions. 

First, I would like to comment on sound-power normalization. As for 
discrete tones, there are numerous mechanisms which one can postulate 
as giving rise to random noise from compressors. While I, myself, and 
others, have pointed to the likely dominances of one or another mechanism 
for a particular class of machine operating in particular circumstances, it 
is unlikely that this will hold universally. Rather there will be transition 
regimes in which a given mechanism gives way to a stronger one. Evidence 
of this is well known in the case of transonic compressor discrete tone noise 
when conditions are reached where rotor-associated force fields may 
propagate and dominate the conventional rotor-stator interaction tones. 
Dr. Morfey has elected to inspect the broadband noise from a wide range 
of compressors (table I of his paper) in terms of frequency scaling using a 
length scale associated with the rotor blade momentum thickness. The 
sound-power level normalizing technique he uses is also heavily de- 
pendent upon the mechanism of lift fluctuations over blade surfaces 
caused by turbulence wakes from upstream rows. 

While this approach might be realistic for multistage compressors (say 
three stages and above), work done over the last 2 to 3 years at Rolls- 
Royce (ref. D-l) indicates that for subsonic single-stage fans, at least, 
other mechanisms can and do strongly influence the radiated sound char- 
acteristics. These include blade interactions with the turbulence structure 
of the intake flow. I would not, therefore, expect Dr. Morfey’s normalizing 
functions to cope with both multi- and single-stage machines, and perhaps 
this is another reason for the tendency of the engine E data to fall at the 
lower scatter band. 

Another subject I wish to discuss is frequency resolution of data. 
The data in figure 2 of reference 5 used by Morfey were frequency-analyzed 
using 100-Hz Spectral Dynamics equipment with approximately 2-second 
averaging (from the 15 ips recorded data on a 30-in. tape loop). 

It is interesting to inspect such data with even greater resolution and 
averaging. An example is the Rolls-Royce Conway compressor noise for a 
well-subsonic operating condition, as shown in figure D-l, where 20-Hz 
Spectral Dynamics equipment was used with 7-second averaging. Besides 
the “sum-and-differencc” tones of the type noted by Crigler, et al. (ref. 
D-2), the signal clearly contains a large range of tones at compressor 
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shaft orders. Most of these protrude less than 5 dB and hence would 
appear much as the data of Snow (ref. 5) when resolved to only 100 Hz. 

The data for engines — A, B, C, and D being of 5-percent bandwidth — 
would also fail to resolve adequately the compressor shaft order signals. 

Although I do not agree that the sound-power level data collapse 
obtained by Morfey is a significant improvement on past attempts, his 
correlation does, however, seem reasonable. Nevertheless I feel that in 
view of my comments above this can only be fortuitous. The directivity 
data collapse, is, however, surprisingly good and I wonder whether Dr. 
Morfey would agree that this could be due to a much closer association 
with discrete tone characteristics than he had assumed. 

Finally, I would like to comment on Dr. Morfey’s optical model for 
high-frequency directivity. A marked difference noted in directivity 
between the single-stage data and the remaining engines is interesting. 
Dr. Morfey has developed a simple optical model to explain this, but I 
wonder if he has attempted to use the approaches of Kaji and Mani 
(refs. 15, D-3). The neglect of flow-associated effects on the propagation 
of sound through blade rows seems unwarranted and I feel that one 
needs an integral approach to the problems of Dr. Morfey’s sections 
entitled “Optical Model for High-Frequency Directivity” and “Effects 
of Source Motion and Mean Flow.” I also feel that the effect of blade 
spanwisc twist cannot be ignored since conditions may not be such as to 
prevent transmission over the entire compressor annulus. 

R. MANI (University of Massachusetts) : The paper is an elegant con- 
tribution to the broadband noise problem in turbomachines. It should 
prove particularly useful as a guide to experimental work. In contrast to 
earlier works on the subject, while the author uses similarity arguments as 
the need arises, there is a welcome emphasis on use of first principles in 
computing the sound power, directivity, energy split, etc. The following 
are some more specific comments: 

(1) In the theoretical development the effect of longitudinal (chord- 
wise) velocity fluctuation is included resulting in a Cl term in/(/? 2 , <*i)- 
In computing the sound-power parameter G, the Cl term was dropped. 
Now in the Smith-House paper, the authors explain the higher normalized 
sound power at part speed as due to the cascade operating at lower 
incidence angle and hence in regimes of higher lift-curve slope. The Smith 
and House argument is based on quasi-steady, real fluid effects and I 
wonder if the retention of the Cl term in /(fh, ai) does not provide a 
more natural explanation, based on unsteady airfoil theory, of the trends 
observed by Smith and House. 

(2) With regard to the optical model for transmission through blade 
rows, the mechanism by which a blade row produces reflected and trans- 
mitted waves is different from the one suggested by figure 6. The blade 
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row behaves like a pair of line diffraction gratings coupled by parallel- 
plate waveguides. The reflected and transmitted waves have orientations 
given by the Tyler-Sofrin formula that the tangential wave number k v 
of these waves differs from that of the incident wave by an integer multiple 
of (2ir/d), being the tangential spacing of the blade row. Thus blade rows 
of completely different stagger angles but of the same tangential spacing 
would give the same orientations of reflected and transmitted waves. 
Actually, at high frequencies the effects are more complex than at low 
frequencies since several higher order transmitted and reflected waves are 
produced. Based on some calculations that I have done on this problem, 
I did conclude that the author’s surmise, to wit that the acoustic energy 
transmission coefficients are low in the incidence regions shown shaded 
in figure G, is probably correct. I believe, however, that the mechanism 
involved here is that the degradation of acoustic energy to wake kinetic 
energy (always present in such a problem when moving medium effects 
arc involved) is high in these incidence regimes at high frequencies. 

(3) To students of the theoretical aspects of the subject, attention is 
specially directed to the appendixes, where considerations involved in the 
transition from the discrete-frequency to the continuous-spectrum case 
are comprehensively dealt with. 

(4) I agree with the author’s approach of gradually building the 
refinements into the analysis. Examples are extension of the Kemp-Sears 
work to include chordwise velocity fluctuations and the effect of the cutoff 
weighting factor on the acoustic efficiency of axial fluctuating forces, etc. 
Perhaps the author will agree that there are further refinements not 
mentioned here which could be; significant. Two examples with regard to 
sound-power estimation are that tangential fluctuating forces would make 
a contribution near cutoff and, secondly, work on estimation of unsteady 
blade forces suggests that Kemp-Sears type estimates are rather in- 
adequate near cutoff. 

MORFEY (author) : The discussers raise a number of interesting 
points and I am grateful for their comments. Since the questions mainly 
fall in the three categories used by Air. House, I have grouped my replies 
under these headings. 

I will first address the subject of sound-power normalization. The fact 
that the single-stage fan (E) falls on the lower edge of the scatter band 
(fig. 3) might well be accounted for by a difference in source mechanisms, 
as Mr. House suggests. Figure 2 shows that, compared with the multi- 
stage machines, the mean spectrum shape for engine FI is deficient at 
high frequencies. 

Dr. Alani’s point about the C L effect in/(ft>,ai) deserves further study. 
Retention of the C L term would raise the normalized level of engine E at 
maximum speed by around 5 dB. Because this machine is a single-stage 
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fan, however, the correction may not be much different over the whole 
speed range. Probably a more important reason for the failure to follow 
a V 6 speed variation (as noted in ref. 3) is the breakdown of Strouhal 
number scaling in the frequency spectrum. 

Dr. Lowson’s final question is also relevant to the normalization pro- 
cedure. A comparison of the dipole and quadrupole mechanisms — with 
a number of simplifying assumptions — has been made for the related 
problem of an isolated rotor in nonuniform flow (sec Paper D3, Lough- 
borough Symposium on Aerodynamic Noise, September 1970). If one 
adopts the Kemp-Sears method (ref. 9) of calculating the unsteady blade 
forces, the comparison shows that dipole noise is likely to predominate 
unless the rotor pressure ratio exceeds 1.1. However, this conclusion may 
have to be modified as theory and experiment yield more data on the 
interaction of blade rows with turbulent flow. 12 

The second subject I will reply to is frequency resolution of the data. 
This is a key question which has not been satisfactorily resolved in 
published data, and the danger of even 100-Hz bandwidth analysis is well 
illustrated by the Conway spectrum. In fact the engine E data used in the 
paper were analyzed initially at 20 Hz, and no shaft order tones were 
detectable at this bandwidth. Nevertheless, the possibility remains that 
the spectra consist of overlapping peaks centered at multiples of shaft 
rotation frequency. These peaks would arise, of course, from irregularities 
in blade geometry or spacing round the rotor disk. 

Mr. House suggests that such “semiordered” noise would show a better 
directivity collapse than true broadband noise; but the evidence from tone 
radiation patterns (ref. 8) indicates that the more ordered the source the 
less correlation there is between engines. 

Finally, I wish to reply to comments made concerning the optical model 
for high-frequency directivity. It would be interesting to take Dr. Mani’s 
theory on sound transmission through blade rows and apply it at very high 
frequencies. I would expect, as he suggests, that the wave theory would 
support the trend predicted by the optical model. If this is the case, it 
might be worth extending the optical model to include the effects of flow 
and blade twist, as indicated by Mr. House. 


12 (Note added in proof) Substantial advances have been made, since this paper was 
presented, in understanding the radiation of sound from blades in turbulent flow. 
Reference D-4 gives an updated review of the situation; it is arguable that the spec- 
trum collapse on L/X should be replaced by a collapse on rotor chord/wavelength. 
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